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Summary. — We review recent results of Lattice QCD calculations relevant for 
flavour physics. We discuss in particular the hadronic parameters entering the am- 
plitudes of K° - K°, D° - D° and B° - B° mixing, the B- and D-meson decay 
constants and the form factors controlling B-meson semileptonic decays. On the 
basis of these lattice results, which are extensively collected in the paper, we also 
derive our averages of the relevant hadronic parameters. 

PACS 12.38.Gc - Lattice QCD calculations. 



The aim of this talk is to illustrate the status of the art of lattice QCD calculations 
relevant for flavour physics and to provide averages of the hadronic quantities which are 
useful for phenomenological analysis of flavour physics, within or beyond the Standard 
Model. These averages will be used in particular by the \JT fit collaboration in the 
unitarity triangle analysis (UTA) [l]-[4]. We believe, however, that the whole discussion 
given in the present contribution can be of more general interest. Even though the 
process of averaging lattice results unavoidably involves some degree of subjectivity, for 
all the quantities discussed in this work we also provide extensive collections of lattice 
results, as well as the references to the original papers. In addition, we try to enrich the 
discussion with some technical considerations about the accuracy of the various lattice 
calculations, which can be useful for the non-expert of the field in order to have a deeper 
view and to draw their conclusions. 

The status of the art of Lattice QCD calculations is represented by unquenched si- 
mulations performed with cither Nf = 2 or Nf = 2 -|- 1 dynamical quarks and values of 
the light (i.e. up/down) quark masses well below ms/^, where rus is the physical strange 
quark mass. Both these features are particularly important to keep the lattice systematic 
uncertainties under control. On the one hand, it is well known that the quenched appro- 
ximation introduces a systematic error which, besides depending on the specific quantity 
under investigation, is also difficult to evaluate (if not by performing the corresponding 
unquenched calculation). On the other hand, lattice simulations performed with suffi- 
ciently light values of the up and down quark masses (ruud < 't^s/2) allow to rely on the 



(*) Speaker 



2 V. LUBICZ and C. TARANTINO 

predictions of chiral perturbation theory when performing the chiral extrapolations of 
the lattice results to the physical masses. 

If the possibility of performing unquenched calculations at light quark masses repre- 
sents a major advance of Lattice QCD in the last few years, it has also to be noted that 
these calculations have not always reached yet the same accuracy in controlling other 
sources of systematic effects as done in the more recent (and less expensive) quenched 
calculations. For instance, some of the unquenched results have not yet been extrapo- 
lated to the continuum limit, and non-perturbative renormalization techniques have not 
always been implemented. Moreover, while most of the quenched results have been de- 
rived by many lattice collaborations, often using different approaches, in several cases 
only few (if any) unquenched results are already available for a given quantity. For in- 
stance, the unquenched results for the kaon B-parameter have been only obtained so far 
at fixed (and rather large) values of the lattice spacing, whereas the matrix elements 
of the full basis of four-fermion operator relevant for K'^ — K^ mixing in generic New 
Physics models have been only obtained so far within the quenched approximation. 

In what follows we illustrate in some detail the results for various lattice parameters 
and present our averages. We will not be able to discuss in the present contribution the 
whole set of hadronic parameters available from the lattice and which are relevant for 
flavour physics. In particular, we are not going to review the accurate lattice results for 
the quark masses, for the ratio of decay constants JkI Jt: and for the vector form factor 
controlling the semileptonic Kii decays. For the latter two quantities, which allow the 
important determination of V^s (the Cabibbo angle), we refer to the review by A. Juttner 
at the Lattice'07 conference [5]. 

1. The kaon i?-parameter Bk 

A collection of quenched and unquenched lattice results for i?^^(2 GeV) [6]- [12] is 
shown in fig. 1 as a function of the square of the lattice spacing a. Among the quenched 
results we have selected all determinations which are C'(a)-improved, i.e. that are affected 
by discretization errors of O(a^) or smaller, and have been obtained at several values 
of the lattice spacing, thus allowing an extrapolation to the physical continuum limit. 
An additional feature of all results shown in fig.l is that they have been obtained with 
lattice discretizations of the fermionic action such that the AS — 2 four-fermion operator 
defining Bk is subject to multiplicative renormalization only. This allows to significantly 
increase the accuracy of the determination of the kaon B-parameter. 

As illustrated in fig.l, the unquenched results for Bk, represented by filled points in 
the plot, have been all obtained so far at a fixed value of the lattice spacing. Thus, at 
variance with the quenched determinations, the unquenched results cannot be extrapo- 
lated to the continuum limit, and the size of discretization errors affecting them cannot 
be easily estimated. On general grounds, we note that there is no reason to expect dis- 
cretization errors in the unquenched case to have smaller size than those affecting the 
quenched results. 

The comparison between quenched and unquenched determinations of Bk obtained 
with the same lattice action and at similar values of the lattice spacing does not show 
any evidence of the quenching effect within the present uncertainties. Moreover, the 
comparison between Nf — and Nf — 2 results obtained with the same lattice action 
but at different values of the lattice spacing, see for instance the case of the DBW2-DWF 
action (empty and filled circles in the plot), suggests that quenched and unquenched 
results could be affected by similar discretization errors. This observation has motivated 
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Fig. 1. - Lattice QCD results for B^ {2 GeV) [6]-[12] plotted as a function of the lattice spacing 
square. All results selected in the plot are 0{a)-improved and required only multiplicative 
renormalization of the relevant AS = 2 four-fermion operator. Filled points are unquenched 
results obtained with either Nf — 2 or Nf = 2 + 1 dynamical quarks. Empty points are quenched 
results. Our average is shown by the horizontal band. 



the RBC collaboration to consider their Nf — and Nf — 2 results together in order 
to perform a combined extrapolation to the continuum limit [7] . By taking into account 
the scale dependence suggested by the quenched results, we quote as final average for 
Bk the value 

(1) B^{2 GeV) = 0.55 ± 0.05 , 

which corresponds to the renormalization group invariant parameter 



(2) 



Bk = 0.75 ± 0.07 



Clearly, a better estimate of discretization errors in the unquenched case will be obtained 
once the results of unquenched simulations performed at different values of the lattice 
spacing, which are currently in progress, will be available. 

Our average in eq. (1) is shown by the horizontal band in fig.l. It can be compared 

with the average Bf^(2 GeV) = 0.58 ±0.03 ±0.06 quoted by Dawson at Lattice'05 [13], 
where the central value and the first error corresponded to the average of the quenched 
results in the continuum limit, while the second error was an estimate of the quenching 
effect. The recent unquenched results for Bk allow to remove the latter uncertainty and 
slightly decrease the previous central value. 



2. The B- mesons decay constants fs^ and fs 

The unquenched lattice QCD results for the i?-mesons decay constants Jb^ and fs 
and for their ratio fssl fs [14]- [20] are shown in fig. 2. They have been obtained by 
treating the h quark on the lattice with two different approaches, either FNAL [21] or 
NRQCD. 
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Fig. 2. - Lattice QCD results for the B-mesons decay constants fss (top left), Jb (top right) and 
the ratio js^l fs (bottom) obtained with Nf — 2 and Nf =2 + 1 unquenched simulations [14]- 
[20]. A star in the legends labels preliminary results. Our averages are shown by the vertical 
bands. 



Our average for the B^ decay constant, 



(3) 



Jb, = 245 ± 25 MeV 



is based on the results shown in the upper pfot of fig. 2. The lattice determination of 
Jb is more delicate, because its value is enhanced by chiral logs effects relevant at low 
quark masses. In order to properly account for these effects, simulations at light values 
of the quark mass (typically mud < Tns/2) are required. For this reason, we derive 
our average for /s and for the ratio fss/fB by taking into account only the results 
HPQCD'05 [19] and FNAL/MILC'07 [20] (see fig.2), which use the MILC gauge field 
configurations generated at light quark masses as low as rUs/S. In this way we obtain 



(4) 



/b.//s = 1.21 ±0.04. 



Finally, we derive directly from the previous averages of /b^ and /bs/ fs the estimate of 
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Fig. 3. - Lattice QCD results for the D-raesons decay constants fos (top), fo (center) and for 
the ratio fo^/fo (bottom) obtained with Nf — 2 and Nf = 2 + 1 unquenched simulations [14, 
16, 20, 23, 24]. A star in the legends labels preliminary results. Our averages are shown by the 
vertical bands. The corresponding experimental averages [22] are also shown for comparison. 



the i?-nieson decay constant 

(5) /b = 200 ± 20 MeV 



3. The Z3-mesons decay constants fn^ and /£> 

The CKM matrix elements Vcs and Vcd, which control the rate of leptonic Df -^ t^vi 
and D+ -^ t^vg decays, are well constrained by the unitarity of the CKM matrix, which 
predicts \Vcs\ — 1 — A^/2 and \Vcd\ — A (up to 0(X'^) corrections), where A is the 
Wolfenstein parameter. By using these constraints, the measurements of the leptonic 
decay rates then allow a determination of the D-mesons decay constants, fos and fo- 
The current experimental averages, as evaluated in ref . [22] , are 



(6) 



/. 



EXP. 



= 273 ± 10 MeV 



/. 



EXP. 



= 205.8 ± 8.9 MeV 



(JdJId) 



EXP. 



1.33 ±0.07 . 



The unquenched lattice QCD results for these constants [14, 16, 20, 23, 24] are shown 
in fig. 3, where they are also compared with the experimental averages given in eq. (6). 

Among the various lattice determinations, which show overall a very good consistency 
among each other, it deserves to be noted the result of the HPQCD collaboration [24], 
which quotes fo^ and fo with the impressive accuracy of about 1% and 2% respectively. 
Remarkably, while the result for /d of ref. [24] is in good agreement with the experimental 
average, the prediction for fjj^ differs from the experimental value by approximately 3 
standard deviations. This discrepancy has been interpreted in [25] as an evidence for 
New Physics effects in leptonic Ds decay. 

A critical review of the lattice study of ref. [24] is not possible yet, since important 
details on the analysis are not given in the paper, being postponed to a longer publi- 
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Fig. 4. - Lattice QCD results for Bg^{mb) and Bgf{mb) (left) and for their ratio BbJEb^ 
(right) [17], [28]-[32]. A star in the legend labels preliminary results. Our averages are shown 
by the vertical bands. 



cation [24]. Of particular relevance, in this respect, are the details of the Bayesian fits 
implemented to perform the combined continuum and chiral extrapolation, as well as 
the knowledge of the precise functional form assumed for this extrapolation. The large 
number of free parameters introduced in the fit (specifically, 45 parameters with 28 data 
points) is the source of some concern, particularly in view of the remarkable accuracy 
claimed on the final results. We also mention that the use of the so called "fourth root 
trick" in dynamical simulations of staggered quarks, which the results of ref. [24] are 
based on, is controversial. This approach has been the subject of intensive theoreti- 
cal and numerical investigation in the last few years. Despite the valuable progresses 
achieved in our understanding of the rooting procedure, the dispute on its validity has 
not been, and will likely never be, completely resolved (for extensive discussions, see e.g. 
the recent reviews [26, 27] at the Lattice conferences). 

Waiting for a confirmation of the precise predictions of ref. [24] by other lattice cal- 
culations, in deriving averages for the D-mesons decay constants we choose to quote 
uncertainties which are larger than those given in ref. [24], and comparable in size to 
those presented by the other results. Thus, we quote as our final averages 

(7) /d, = 250 ± 15 MeV , /d = 212 ± 14 MeV , /^J/d = 1.18 ± 0.03 , 

which are the values also shown by the vertical bands in fig. 3. 



Matrix elements for -Bj]/^ 



-Bj]/^ mixing 



A collection of quenched and unquenched lattice results for the Bd and Bg bag pa- 
rameters and for their ratio [17], [28]- [32] is shown in fig. 4. We firstly observe that the 
dependence on the light quark mass, that should allow to distinguish between Bj, and 
J5s, is practically invisible. Moreover, the unquenched results, obtained with Nf = 2 and 
Nf = 2 + 1 dynamical quarks, tends to be slightly lower than the quenched determina- 
tions, though still well compatible within the errors. On the basis of the unquenched 



FLAVOUR PHYSICS AND LATTICE QCD: AVERAGES OF LATTICE INPUTS FOR THE UTA 7 

results we quote the average 

(8) B^{mb) - B^inib) = 0.80 ± 0.08 , 

which is illustrated by the vertical band in fig. 4 (left) and corresponds to the rcnormali- 
zation group invariant parameters 

(9) Bb^=Bb,, = 1.22 ±0.12 . 

For the ratio Bg^/Bg^, where statistical and systematic uncertainties partially cancel 
(see fig. 4 right), we quote 

(10) BbJBb^ ^ i.oo±om . 

Combining the above results with the averages quoted for the _B-mesons decay con- 
stants, eqs. (3)-(5), we then obtain 



/s. V Bb^ = 270 ± 30 MeV , /b V ^s, = 225 ± 25 MeV , 



(11) C= ^^^= = 1.21 ± 0.04 



fB\^BB, 



It is important to compare the average of Jb^v Bb^ of cq. (11), obtained from Jb 



and Bb, separately, with the direct calculation of Jb^ y Bb^ from the whole i?,5-mixing 
matrix element. For the latter, the unquenched results obtained in JLQCD'03 [17] and 

HPQCD'06 [31] are /b,^/^ = 245 ± lOti? MeV and /s.^^ = 281 ± 21 MeV 
respectively, well consistent with the average of eq. (11). 

5. Vcb exclusive 

The lattice results for the zero momentum transfer form factors of -B ^ 0*1^ and 
B -^ Dlv semileptonic decays [33]-[38], denoted as F{1) and G{1) respectively, are 
shown in fig. 5. In the infinite mass limit of both the charm and bottom quarks these 
form factors are normalized to unity, up to small radiative corrections. Fig. 5 shows that 
the percent level accuracy required to evaluate the l/mg corrections to the static limit 
has been reached by the lattice calculations. Notice that the form factors for both the 
B -^ D*lv and B —=> Dlv decays have been determined by the FNAL and ROMA- 
TOV collaborations using different techniques, and the agreement between the results 
is reassuring. From a phenomenological point of view, the semileptonic decay in the 
vector channel, B -^ D*lv, plays a privileged role, since the experimental accuracy in 
the determination of the decay rate is higher by almost a factor 10 than the one reached 
for B -^ Dlv decays. 

For the form factor F{1), an unquenched determination with Nf = 2 + 1 dynamical 
quarks has been recently obtained by the FNAL group, and it is in good agreement with 
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Fig. 5. - Lattice QCD results for the form factors F{1) (left) and G(l) (right) [33]-[38] controlhng 
B -^ D*lu and B — > Dlv decays. A star in the legend labels preliminary results. Our averages 
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the quenched estimates by the FNAL and ROMA-TOV collaborations (see fig. 5 right). 
Thus, we quote the unquenched result as the best estimate of the form factor, namely 



(12) 



F(l) = 0.924 ±0.022 



In the case of i? ^ Dli/ decays, the unquenched result for G{1) by the FNAL collabo- 
ration has been only presented at the Lattice 2004 conference [38] and it has not been 
submitted for publication. Therefore, in this case, we prefer to quote as our final estimate 
the value 



(13) 



G(l) = 1.060 ±0.035 



obtained by combining the unquenched and the two quenched determinations by the 
FNAL and ROMA-TOV collaborations, also shown in fig.5 (left). 

In order to extract Vcb, we use the HFAG experimental averages \Vcb\ ■-F'(l) = (36.18± 
0.55) • 10-3 and \Vcb\ ■ G{1) = (42.3 ± 4.5) • lO^^ [39], thus obtaining 



(14) 



\Vcb\ (excl. 
\Vcb\ (excl. 



(39.2 ±1.1)- 10" 
(39.9 ±4.4) -10" 



from B 
from B 



D*lv 
Dlv . 



Clearly, the determination of Vcb is dominated by the result obtained from B 
decays, which we thus take as the best estimate for the CKM matrix element. 



D*lv 



— Vub exclusive 



Lattice QCD results for the form factor controlling the exclusive semileptonic B -^ irli/ 
decay are usually considered to evaluate the integrated decay rate in the large momentum 
transfer region, which is the one directly accessible to lattice calculations. In fig. 6 we 
present a compilation of both quenched and unquenched lattice results [38], [40]- [44] for 



the quantity FF{q^ > 16 GeV^) = r{q'^ > 16 GeV^)/|Kb| ■ They have been obtained 
using different approaches to treat the b quark (relativistic QCD, NRQCD and FNAL). 
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A very good agreement among all the results is observed, and the effect of the quenching 
error is not visible so far within the present uncertainties. 
By averaging the lattice results shown in fig. 6 we obtain 



(15) 



FFiq"^ > 16 GeV^) = 2.04 ± 0.40 ps" 



which can be combined with the experimental branching fraction BF{q^ > 16 GeV ) = 
(0.38 ± 0.03 ± 0.03) • 10-4 quoted for B° by HFAG [39] and with the lifetime Tgo = 
1.530 ± 0.009 ps to obtain 



(16) 



iKbl (excl.)^Qc^ = (35.4±4.0).10- 



For these decays also the information coming from QCD sum rules on the small 
momentum transfer region can be considered. Ref. [45] provides FF{q^ < 16 GeV^) — 
5.44 ± 1.43 ps"^ which, combined with the HFAG average, BF{q'^ < 16 GeV^) = (0.95 ± 
0.05 ±0.05) ■ 10""' and the S" lifetime quoted above, gives |y„b| {Gxc\.)QQ]j_gj^ = (34.7 ± 
4.8) • 10~''. The more recent and updated analysis of ref. [46] quotes 



(17) 



\Vub\ (excl.) 



QCD-SR 



(35±4±2±1) -lO^^ , 



in good agreement with the previous QCD sum rules result and the lattice determination. 
We thus quote as our final average 



(18) 



\Vub\ (excl.) = (35.0 ±4.0) -10- 



7. Matrix elements for K'^ — K^ mixing: full basis of four-fermion operators 

In models of physics beyond the Standard Model, the effective Hamiltonian which 
describes the K^ — K^ mixing amplitude involves in general the complete basis of AS" — 2 
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four-fcrmion operators. A common choice for this basis is constituted by the operators 

Oi - 5"74l - 75)rf"5^7M(l - 75)^'' , 

(19) 03 = 5"(l-75)d^5^(l-75)rf", 

04 = 5"(l-75)d"s'3(l+75)ci'^, 

05 = 5"(l-75)d^5^(l + 75)d", 

where a and /3 are colour indices, together with the operators Oi,2,3 obtained from 
Oi,2,3 with the exchange 75 -^ —75- In chirally invariant renormalization schemes, the 
operators Oi have the same matrix elements of the Oi and, for this reason, they will not 
be further discussed in what follows. 

Omitting terms which are of higher order in chiral perturbation theory, the B- 
parameters are introduced using the expressions 

{K"\Orif^)\K") = ^My],Bfi^,) , 

/ \ 2 



(20) (A-"|0,(M)|A-") ^ i ( „..(,'^';„,(,) )'"i-/i-«f(rt . 

(A-°|0.(m)|A'°) = 2 (" "*" , . )%4/|Bf(rt . 

3 \ms(ii) +md[p.) J 

where Sf'^(^) == B/^ (a*). 

The matrix elements of the full operator basis for K'^ — K'^ mixing have been com- 
puted, at present, only in three lattice studies [47, 48, 49], all performed within the 
quenched approximation. Refs. [47] and [48] use tree-level improved Clover and overlap 
fermions respectively, and both implement non-perturbative renormalization with the 
RI-MOM method. The results of ref. [49] are obtained by using domain wall fermions 
and renormalizing the operators with one-loop perturbation theory. They have been pre- 
sented at the Lattice 2006 conference but a final analysis has not been published. The 
values of the _B-parameters from refs. [47]-[49] are collected in table I (top) (^). In table I 
(bottom) we collect, instead, the corresponding values of the ratios of the non Standard 
Model to the Standard Model matrix elements, 

.rf. ^ {K^\Ot{^i)\K^) 



( ) Wc have selected from ref. [47] the results obtained at the largest values of the lattice 
coupling, and from ref. [49] those obtained with the Iwasaki gauge action on the larger volume. 
In the latter case, we add to the statistical error the uncertainty coming from the difference 
between the results obtained with the Iwasaki and the Wilson plaquette gauge actions, which 
provides an estimate of O(a^) discretization effects. 
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Bf 


Bf 


Bf 


Bf 


jDsd 

B5 


Ref. [47] 


0.68(21) 


0.67(7) 


0.95(15) 


1.00(9) 


0.66(11) 


Ref. [48] 


0.56(6) 


0.87(8) 


1.41(16) 


0.94(6) 


0.62(8) 


Ref. [49] 


0.52(4) 


0.54(2) 


0.71(2) 


0.70(1) 


0.62(1) 














Rf 


Rf 


Rf 


Rf 


Rf 


Ref. [47] 


1 


-7(2) 


1.9(6) 


12(3) 


2.6(9) 


Ref. [48] 


1 


-16(3) 


5.2(9) 


21(3) 


4.6(9) 


Ref. [49] 


1 


-19(1) 


5.0(3) 


30(3) 


8.8(7) 



Table I. - B-pararaeters (top) and ratios Rf' (bottom) for the full basis of four-fermions 
operators in K^' — K'^ mixing. Results are given in the RI-MOM scheme at the scale fi = 2 GeV. 



From the results shown in table I we observe that the differences among the three 
determinations of refs. [47]- [49] are typically larger than the quoted uncertainties. Clearly, 
new lattice studies of these matrix elements would be necessary in order to clarify the 
situation. The observation that the discrepancies are much more pronounced in the case 
of the ratios Ri might suggest that some of the systematic errors cancel, at least in part, 
in the determination of the S-parameters. For this reason, we average the results for the 
B-parameters and quote: 



(22) 



Bf = 0.7(2) , Bf = 1.0(4) , Bf = 0.9(2) , Bf = 0.6(1) 



in the RI-MOM scheme at the scale /i = 2 GcV. For the parameter Bf^ = Bk we refer 
instead to the average given in cq. (1) which, translated into the RI-MOM scheme at 
PL = 2 GcV, corresponds to Bf = Bk ^ 0.54(5). 

In order to derive from the average values of the B-parameters the ratios Rf^, a 
determination of the strange and down quark masses is required. We use, at fi ~ 2 GeV, 



RI 



,RI 



(m 

way we find 



) — 135 ± 18 MeV, corresponding to (m 



MS 



,MSn 



110±15MeV. In this 



(23) Rf = -11±4 



i?f = 3.1 ±1.5 



Rf = 17 ± 6 



Rf = 3. 



±1.3 



8. 



Matrix elements for B'^,^- 



-B^,^ mixing: full basis of four-fermion operators 



In the case of B mesons, besides the matrix element of the operator Oi which enters 
the mixing amplitude, also the operator O2 (sometimes denoted as Os) is of interest 
in the Standard Model, since it contributes to the theoretical prediction of the lifetime 
difference AF of the neutral B mesons. As in the case of K^ — K^ mixing, in generic 
New Physics models the effective Hamiltonian describing /S.B = 2 transitions receives 
contribution from the full basis of four fermions operators. Here we define the operators 



12 



V. LUBICZ and C. TARANTINO 





B- 


B2' 


i3- 


Bf 


Bt' 


Ref. [29] 


0.88(5) 


0.84(4) 


0.91(9) 


1.15(6) 


1.74(7) 


Ref. [30] 


0.86(5) 


0.86(5) 


— 


— 


— 


Ref. [31] 


0.76(11) 


0.84(13) 


0.90(14) 


— 


— 



Table II. - B-parameters for the full basis of four-fermions operators in B^ — 5° mixing. The 
resuhs are given at the scale /i — nit in the MS scheme of ref. [50] for B1-B3 and of ref . [51] for 
B4-B5. 



O1-O5 for B' 



d/s 



^d/s riiixing and the corresponding B-paranieters as in eqs. (19) and 
(20), with the obvious replacements s —^ b and d -^ d/s. 

Lattice results for the full basis of AB = 2 four-fcrmion operators have been obtained 
only in ref. [29] , within the quenched approximation. A quenched determination of Bi 
and B2 has been also performed in ref. [30] (^). Recently, the matrix elements of O1-O3 
have been also computed in a simulation with Nf = 2 + 1 dynamical fermions [31]. All 
these results are collected for the Bs system in table II, with the i?-parameters defined 
at the scale fx — 2 GeV in the MS scheme of ref. [50] for B1-B3 and of ref. [51] for 
B4-B5 (these are the renormalization schemes usually considered for these quantities). 
Note that the larger errors quoted by ref. [31] on the B-parameters, with respect to 
refs. [29, 30], are due to the fact that the i3-parameters are extracted in [31] from the 
values of the corresponding matrix elements, rather than directly calculated from the 
more precise ratios of 2- and 3-point correlation functions. On the other hand, in order 
to evaluate from the B-parameters the full matrix elements of four-fcrmion operators, 
which are eventually the quantities of physical interest, the results of refs. [29, 30] for 
the B-parameters have to be combined with the value of the decay constant fss , which 
thus increases the final uncertainty. 

The SU{3) breaking effects for the full basis of B-parameters have been computed in 
ref. [29], that quotes 



(24) 



{B^/B 



bd\ 



,5 = {0.99(2), 1.01(2), 1.01(3), 1.01(2), 1.01(3)} 



These ratios are well compatible with unity within the errors. Thus, we choose to quote 
common averages for Bd and Bs mixing. 

The results collected in table II show a very good agreement among the central values 
of the three (two) determinations of B2 (B3). On the other hand, the comparison of 
several lattice determinations for the B-parameter Bi = Bb has led us to quote in 
eq. (8) a final uncertainty on this quantity of about 10%. Therefore, we choose to add 
the same 10% of systematic uncertainty also to the final estimates of B2 and B3, as well 
as to those of B4 and B5 for which only a single (and rather old) lattice calculation exists. 



(^) Ref. [30] quotes final results for the B2 parameter only, although both the matrix elements 
of the operators O2 and O3 have been evaluated, since they mix under renormalization. 
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In this way, we quote 

(25) bI'^ = 0.85(10) , bI'^ = 0.90(13) , B^^ = 1.15(13) , Bg"^ = 1.74(19) , 

for both q = d,s. The results in eq. (25) are given at the scale jjl = rrii, in the MS scheme 
of ref. [50] for B2-B3 and of ref. [51] for B4-B5. For the parameter B^"^ = Bb^ we refer 
instead to the average given in eq. (8) (in the same MS scheme of ref. [50]). 

The values of the i3-parameters in eq. (25), combined with the average of the Bs- 
meson decay constant of eq. (3), allow to evaluate the quantities 



(26) 7^,(mb) = ( _ ,_ r,""! ^^) ) fBs ^B\^{nn,) , z - 2, 



Me 
mb{mb) +ms{mb) 



The squares Tif are proportional to the corresponding operator matrix elements. Us- 
ing m^^{mb) = 4.21(8) GeV and mf^{mb) = 87(12) McV (which corresponds to 
mp(2 GeV) = 105(15) MeV), we obtain 

712 = 282±34 McV , 7^3 = 290±37 MeV , 7^4 = 328±39 MeV , 7^r, ^ 404±47 McV . 

(27) 

For 7^2 ^nd TZ^ these results are well consistent with those quoted by the unquenched 

calculation of ref. [31], namely 7^2 ^ 295 ± 22 MeV and 7^3 == 305 ± 23 MeV. 

9. Matrix elements for I?'^ D'^ mixing 

The matrix elements of the full basis of four-fermion operators for D'^ — D^ mixing can 
be obtained as a byproduct from ref. [29]. This work provides quenched results for the 
B-parameters corresponding to heavy-light meson masses equal to Mp = 1.75(9) GeV 
and Mp = 2.02(10) GeV respectively. By interpolating between the two sets of results, 
one obtains for the physical D-mesons the values i?™|i=i,...,5 — {0.85(2), 0.82(3), 1.07(5), 
1.10(2), 1.37(3)}, in the RI-MOM scheme at the scale ^ = 2.8 GeV. The S-parameter 
Bo = B1" has been also computed in ref. [52], always within the quenched approximation 
and using non-perturbative RI-MOM rcnormalization, with domain wall fermions. The 
result is affected by a large uncertainty: they obtain _B™ = 0.85(2)(11) in the MS scheme 
at /^ = 2.0 GeV. Translated into the RI-MOM scheme at the scale ^ = 2.8 GeV, this 
value corresponds to _B™ = 0.81 (2) (10), well consistent with the result of ref. [29]. 

By considering that the results of ref. [29] are obtained from a single quenched lattice 
calculation, we add in the final averages a systematic uncertainty of 10%, as we have 
done in the case of i?° — i?° mixing. In this way we obtain 

B™ = 0.85(9) , S^" == 0.82(9) , B^" = 1.07(12) , 
(28) B™ = 1.10(11) , B™ = 1.37(14) , 

in the RI-MOM scheme at the scale ^i == 2.8 GeV. 
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